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ABSTRACT 


This paper documents an invited presentation given to The Boeing Company, Seattle, 
Washington on September 9, 1997. The audience consisted of structural dynamic and flight test 
engineers from the Boeing Commercial Airplane Group who were interested in discussing 
research which may be applied to future flight flutter test programs. A method to compute robust 
flutter margins is described which is a significant departure from traditional methods. This 
method uses the structured singular value, (I, to compute a flutter margin which directly accounts 
for modeling errors such that a worst-case flutter margin is computed with respect to those 
errors. This method may be applied in several ways. A post-flight application uses data sets from 
multiple test points to compute worst-case flutter margins and a worst-case flight envelope. An 
on-line implementation computes flutter margins at each test point to track the flutter margins 
during a flight test. This on-line implementation is the basis for a flutterometer flight test tool 
that displays the distance to flutter at a given test point. Such a tool was not previously possible 
using traditional flutter flight test analysis methods. The F/A-18 System Research Aircraft was 
used to demonstrate these applications using flight data recorded from test points throughout the 
flight envelope. 


NOMENCLATURE 


C damping matrix 

K stiffness matrix 

M mass matrix 

P plant transfer function 

p set of plant transfer functions 

q dynamic pressure 

Q unsteady aerodynamic force matrix 

s Laplace variable 

P pole in Pade approximation to a lag 

8- perturbation to dynamic pressure 

A uncertainty operators 

A set of uncertainty operators 

Ap uncertainty in P 

r( state vector 

p structured singular value 

& set of stable, linear time-invariant transfer functions 

oo 

set of finite square-integrable measurements 
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(i. Method : Research Team at Dryden 


- Marty Brenner 

• Initiated research program 

• Main partner on applying robustness to aeroelasticity 

- Len Voelker 

• Formulated initial concept for flutterometer 

• Computed p-k flutter analysis to compare with [i method 

• No living human has more practical knowledge of flutter 

- Larry Freudinger 

• Developed on-line implementation concepts 

- DaveVoracek 

• Chief engineer on F/A-18 SRA who generated and analyzed data 

- Roger Truax 

• Generated finite element model of F/A-18 SRA 

- Mike Kehoe 

• Supervised team efforts and considered practicality issues 


The structural dynamics group at NASA Dryden Flight Research Center 
has been actively involved in this research project for some years. 

The diversity of the team enables significant research issues to 
be addressed by engineers with specialization in that area. 
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(I Method : Publications 


- R. Lind and M. Brenner, Robust Flutter Margin Analysis that 
Incorporates Flight Data , NASA Technical Paper, in preparation. 

- R. Lind and M. Brenner, “Robust Rutter Margins of an F/A-18 Aircraft 
from Aeroelastic Right Data”, AIAA Journal of Guidance, Control and 
Dynamics , Vol. 20, No. 3, May-June 1997, pp. 597-604. 

- R. Lind and M. Brenner, “Utilizing Flight Data to Update Aeroelastic 
Stability Estimates”, 1997 AIAA Atmospheric Flight Mechanics 
Conference , New Orleans LA, AIAA-97-3714, August 1997. 

- R. Lind and M. Brenner, “A Worst-Case Approach for On-Line Rutter 
Prediction” International Forum on Aeroelasticity and Structural 
Dynamics , Rome Italy, June 1997, Vol. 2, pp. 79-86. 

- R. Lind, M. Brenner and L. Freudinger, “Improved Right Test 
Procedures for Flutter Clearance” International Forum on Aeroelasticity 
and Structural Dynamics , Rome Italy June 1997, Vol. 3, pp. 291-298. 

- R. Lind and M. Brenner, “Worst-Case Flutter Margins from F/A-18 
Aircraft Aeroelastic Data” 1997 AIAA Structures, Structural Dy namics 
and Materials Conference , Orlando FL, AIAA-97-1266, pp. 738-748. 


There are several refereed publications which discuss aspects of 
the method for robust flutter margin analysis. The NASA Technical 
Paper, which is expected to appear in early 1998, is the most 
detailed and complete reference document. 







Outline 


• Motivation 

• [i Method for Robust Flutter Margin Analysis 

- Concept of robustness and uncertainty for flutter margins 

- Incorporating flight data 

• Worst-Case Flutter Margins : Post-Flight Analysis 

- Implementation and approaches to utilize flight data 

- Application : F/A-18 SRA 

• Flutterometer : On-Line Analysis 

- Implementation and approaches to utilize flight data 

- Application : F/A-18 SRA 

• Conclusions 

- Research extensions 

- Open discussions 


This presentation has essentially three large sections. Some initial 
comments are given to briefly discuss flutter analysis issues and 
provide the motivation for this research. The main section 
introduces the method for flutter analysis, called the (1 method, 
which produces robust flutter margins. This portion of the 
presentation details the issues of robustness and demonstrates 
their applicability to flutter analysis. The other large sections 
describe implementations of the fi method and demonstrate 
these implementations on the F/A-18 Systems Research Aircraft. 
The first implementation is a post-flight analysis which computes 
worst-case flutter margins to define a flight envelope. The second 
implementation is the on-line tracking of robust flutter margins 
via the flutterometer concept. 
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Motivation 


Consider current flight flutter test procedures 



A test point 


decision 


data 






damping 

estimate 


The procedures for flight flutter testing to clear a flight envelope are 
fairly standard throughout NASA and the flight test community, 
although implementation details vary widely between organizations. 
Response data is recorded from the aircraft at a stabilized test point 
and telemetered to the control room. The data is analyzed using 
strip chart monitoring and several computational algorithms in 
both the time and frequency domains. Estimates of the modal 
damping ratios and trends for those ratios as the envelope is expanded 
are used to determine the next suitable test point which does not incur 
an excessive level of risk. 
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Motivation 


NASA Dryden research to improve flight test efficiency 
by utilizing flight data 



confident 

flutter margin 


NASA Dryden Flight Research Center is researching methods to 
improve flight test efficiency. Increasing efficiency implies 
reducing time and cost along with maintaining a high level of safety 
for the aircraft and crew. The entire scope of flight flutter testing 
is being investigated. The first step in the process investigates better 
signals and mechanisms to excite the aircraft and improve the dynamics 
observed with response data. An RBNB concept is being implemented 
to distribute this data to a variety of analysis facilities. Improved 
signal processing algorithms such as wavelet filtering are being 
utilized to accurately process transient response data. The last step 
in the process is to analyze this data to produce a confident 
flutter margin rather than the damping estimate which is currently 
computed. This last step is the focus of this presentation and 
is addressed by a parameter called p. 
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Motivation 


• Introduce concept of a flutterometer 

- flight test tool 

- indicates distance to flutter 


-10 kft 



i i warning 
■■■ danger 


• Information about flight condition at flutter 

- altitude 

- dynamic pressure 

- airspeed 


These research areas are all steps towards the concept of a 
flutterometer which was envisioned in the 1980's at NASA Dryden 
Flight Research Center. This flutterometer is a flight test tool 
that indicates some measure of distance to a flutter condition. 

The center box in the dial present the altitude at which flutter occurs. 
The values along the dial present the difference between that altitude 
and the current altitude at which the aircraft is flying. 

The type of measure, such as altitude or dynamic pressure or airspeed, 
can be selected to match units desired by the pilot or engineer. 

The main point is this tool provides a quantitative value of the flight 
conditions at which flutter occurs. This tool can drastically 
increase flight test efficiency since test points can be safely chosen 
with greater rate of expansion of the flight envelope. 



Motivation 


M*$A 


Analytical Predictions 

- p-k from FEM and CFD 

- 1 st order perturbations 

- stochastic robustness 



On-Line Estimates 

- damping 

- parameter identification 

- modal filter 



well-developed model — STRENGTHS— ► data describes true aircraft 


model is approximation ^-DRAWBACKS -► can’t extrapolate to margin 


This flutterometer concept can not be effectively implemented using 
any current flutter analysis method. One explanation for this 
is seen by dividing all current methods into two basic categories 
called Analytical Prediction methods and On-Line Estimation methods. 

The Analytical Prediction methods, of which the p-k method is 
the most common, utilize a computational model with no direct consideration 
of flight data. The On-Line Estimation methods, of which tracking 
damping estimates is the most common, utilize the flight data alone. 

Each of these methods has strengths and weaknesses that are directly 
related. The drawbacks to the Analytical Prediction methods are 
eliminated by the On-Line Estimate methods; however, the On-Line Estimate 
methods introduce their own drawbacks which are not problematic 
for Analytical Prediction methods. 



Illll Limitations in Proposed Approaches 

• Analytical Prediction Methods 


- 1st order perturbation analysis 

• perturbation structure may not be realistic 

• flutter margins may be overly conservative 

(Becus,Poiron) 

- Stochastic robustness 

• expensive Monte Carlo simulations 

(Stengel) 

• robustness levels are statistical with no guarantees 

• On-Line Estimation Methods 

- Parameter estimation 

• computationally expensive 

• no convergence or optimality guarantees 

• poor performance for low SNR data 

(Nissim,Feron) 

- Modal filtering 

• model based method is not adaptive or robust 

• possible problems for dense modal spaces 

• its a filtering, not a processing, algorithm 

(Shelley ,AHemang) 


The most common methods of flutter analysis, namely the analytical 
p-k method and the on-line estimation of damping, have been 
recognized as deficient for many years. Several new methods are being 
investigated by various researchers to replace these traditional methods. 
These new methods still fall into the two basic categories and thus 
have the same strengths and drawbacks. These are additional 
limitations that should also be considered for these new methods. 
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I x Method : Information Flowchart 


• NASA method combines strengths of each method 

- Analytical model provides nominal model 

- Flight data provides errors in model 



| P-k 


CD 

Prony | 


Analytical 

Prediction 




Flight Test 
Prediction 


Flight Test 
Estimate 


NASA Dry den Flight Research Center is developing a novel method 
for flutter analysis that combines the strengths of both categories of 
traditional flutter analysis. This method is essentially model based 
so it has the desired predictive nature of an Analytical Prediction 
method; however, this method, unlike traditional Analytical 
Prediction methods, can also utilize the flight data to obtain the 
desired accuracy of the On-Line Estimate methods. This method, 
referred to as the (J. method, introduces a new category of analysis 
called Flight Test Prediction methods. 
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jll Method : Robust Stability 


• jj. uses mathematical representations 

- operators : stable linear time-invariant transfer functions 

- signals : finite bounded square- integrable measurements 


• ji uses Linear Fractional Transformation framework 

- series and feedback interconnections of operators 

- multiple LFT’s with unstructured operators 

result in single LFT with structured operators 


P={f,(P, A) : | A|| < l} 



The (l method is able to make strong claims about robustness 
by utilizing a well developed mathematical framework. The 
underlying concepts are derived from functional analysis and 
introduce abstractions such as ^ operators and signals. 
These abstractions are readily interpreted as standard systems 
concepts such as transfer functions and data measurements. 
The Linear Fractional Transformation is used to represent 
series and feedback operators in a single unified framework. 
The concept of a set of plants, denoted F U (P,A), will be used 
extensively throughout this presentation and refers to the block 
diagram in the Figure. 
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jn Method : Robust Stability 


• Model should be robustly stable to uncertainty A 

- Parametric uncertainty describes errors in specific elements 
• Consider plant with uncertain pole 

? 4»'bU,) : ' s[231 } 


Express this plant in LFT form with A 
P = \f„(P, A): P = 


-L 

irn 

i 

(i+1)<j+2.5) 


j+2.5 

] 


(j+lX$+2.5) 


L1a|<i,A€^ 


- Dynamic uncertainty is a more general type of variation 

• Multiplicative Uncertainty p = {p(7 + A): |A|~ < l} 

• Additive Uncertainty P = {P + A : ||A|— < l} 


The concept of robustness is frequently used in engineering 
terminology to loosely refer to stability and performance of a system 
despite some concept of perturbations. The |i framework 
utilizes a formal definition of robustness and associated perturbations. 
The perturbations are represented by a set of norm bounded operators, 
A, which affect the plant, P, through a feedback relationship 
using the LFT framework. This uncertainty can be parametric 
uncertainty which affect specific elements of a system or 
general dynamic uncertainty which affects groups of signals. 



(i Method : Robust Stability 


• What does robust stability of P mean? 

- No perturbation A with llAlkl can destabilize P 

- Stability margin of P is greater than the size of A 

- Conservative condition for robust stability is llPlI < 1 

• Define the structured singular value jx 

l 

4(P)_ min {a(A) : det(I-PA) = 0} 
AgA 

- multivariable gain/phase margin 

- relates smallest destabilizing perturbation 



• ji is an exact measure of robust stability 

- uncertainty is weighted so (1=1 is the desired condition 

- upper bound computed via convex optimization 


The formalized concept of robustness ensures the plant will be 

stable for any perturbation operator A contained within 

the set A. This concept is directly related to the concept 

of a stability margin. The structured singular value, (X, 

is defined as a necessary and sufficient condition to exactly 

compute the robustness of a plant operator. |i can be interpreted 

as a multivariable gain and phase margin but perhaps the most 

straightforward interpretation is that of the smallest destabilizing 

perturbation. The plant is usually weighted such that the set of 

uncertainty operators has a unity norm bound and the desired robustness 

condition is fi=l. 
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fj, Method : Robust Stability 


• Model validation considers data and model 

- Compares expected measurements with actual measurements 

- Provides measure of plant accuracy 

- Actually tests if model is not invalidated 

since model can not be truly validated by a finite amount of data 


Question : Is there a A in the set A such that F U (P,A) 

could generate measured y in response to w 


• |i can be used as a model validation criterion 

- Form matrix P = f (P,u t y,a>) 

- The model is not invalidated if (a(P)> 1 


A difficulty in using |i to compute robust stability lies in 
choosing a set of meaningful uncertainty operators. Model validation 
algorithms are formulated to address this issue by comparing 
measured data to the set of plants generated by F U (P,A). 

The model validation criterion determines there exists some 
A in the set A such that the plant F U (P,A) 

could have generated the measured data. This criterion actually 
ensures the model is not invalidated since a model can never by 
truly validated by a finite set of data; however, even this restricted 
condition is useful to ensure the uncertainty set is reasonable. 
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j i Method : Robust Stability 


mm-: 


WAS A 




• [i model validation is transfer function based 

- experimental data must fall between upper and lower bounds 

- frequency-varying bounds are determined by uncertainty 

- damping is not estimated 



• is used to consider all A blocks simultaneously 


A approach to model validation in the p framework is 
developed using frequency domain transfer function data. The basic 
concept considers upper and lower bounds for the set of plant models 
with uncertainty and determines if the measured transfer function 
lies within those bounds. This approach uses a simple test at 
each frequency point without requiring estimates of specific system 
and modal parameters. 
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fi Method : Flutter Applications 

W r 


• p has well-known interpretations for controls 

- Gain and phase margins 

- Closed-loop performance despite actuator and sensor noise 

• p has interpretations also for aeroelasticity 

- Describes traditional flutter margin similar to p-k analysis 

- Describes robust flutter margin for plants with modeling errors 

• p Method is valuable tool for flutter analysis 

- Method incorporates flight data 

- Margins are worst-case with respect to some uncertainty 

- p replaces poorly behaved damping as a stability margin 


p was originally developed for control design and analysis and is a well-known 
concept in that field of engineering. NASA Dryden Flight Research Center is 
developing applications of p for aeroelasticity research and has found several 
interpretations for this field, p can be formulated to compute a flutter margin 
which is similar to that generated by traditional analysis methods such as the 
p-k method, p can also represent a robust flutter margin that considers 
the worst-case flutter margin with respect to modeling errors and perturbations 
represented by the set A. 

p can be especially valuable for flutter analysis. The confidence in the 
computed flutter margins are high due to the inclusion of flight data which 
describes the true aircraft. The formalized derivation of p lends some 
mathematical guarantees that the computed robust margins are worst-case 
with respect to a set of perturbations. Also, p is continuous and smooth 
with flight condition variation so replaces the poorly behaved nonlinear 
damping trend parameter to track during a flight test. 
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p Method : Nominal Flutter Margins 


• Consider computing traditional flutter margins 

- traditional margins only use the theoretical plant 

- traditional margins do not consider uncertainty (nominal) 


Question : How can |i represent a traditional flutter margin? 

- (i. relates destabilizing perturbation 

- Flutter margin relates unstable flight condition 


Answer : (i is destabilizing perturbation to flight condition 


The first step in developing the |i method is to generate 
flutter margins which are equivalent to traditional margins 
generated by the p-k method. These margins do not account 
for any uncertainty and are denoted as nominal margins in the 
vernacular. |i is formulated as a flutter margin by 
considering the smallest destabilizing perturbation to a flight 
condition. 
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|i Method : Nominal Flutter Margins 


Consider equation of motion for state vector rj 
Mfj + Cr\ + Kri + qQ(s)rj= 0 


Structural 

Dynamics 


{ M - mass matrix 
C - damping matrix 
K * stiffness matrix 


Unsteady 

Aerodynamics 


{ 


Q - forces matrix 
q - dynamic pressure 


Represent unsteady aerodynamics as transfer function 
a Q b Q I 

Cq Dq J = Dq + Cq(sI-Aq) ! Bq 



The equations of motions for an aeroelastic system utilize a 
structural model and an unsteady aerodynamic force model. The 
unsteady aerodynamic forces can be represented as a finite-dimensional 
state-space model using several standard algorithms. 
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(I Method : Nominal Flutter Margins 


Consider parametric perturbation in dynamic pressure 

q - qo+bq 


Separate perturbation from nominal dynamics 
0 = Mr\ + Cf\ + Kr\ + qQy\ 

= Mr\ + Cf) + Kt | + {q o + + Dq n) 

= Mr\ + Cr\ + At| + q 0 (CQX + Dq r | )+ &i(Cqx + Dgf|) 
= Mr\ + Cr\ + Kr\ + q 0 \CQX + A/r|)+ 

= M^ + Ci] + Kr\ + qo{CQX + Dor\)+w 


Formulate LFT for (i analysis as F u (P 0 ,8 q ) 


0 

1 

0 

0 


-M~ l (K+qoDQ) 

-iVr‘c 

-KC'qoC.Q 



Bq 

0 

Aq 

0 

x_ 

L ^ 

0 

Cq 

0 

w 



The dynamic pressure effects the equations of motion in a linear 
manner so perturbations to this flight condition parameter may 
be easily represented with a feedback operator. The basic 
procedure is to isolate the perturbation from the known nominal 
dynamics and replace this perturbation with a norm bounded operator. 
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|i Method : Nominal Flutter Margins 


• |i directly computes nominal flutter margins 

- range of dynamic pressures is treated as uncertainty 8 q 

- )i considers stability over all uncertainty (dynamic pressures) 


Nominal Flutter Question : 

- What is the smallest q for which P(q) is unstable? {p-k ) 

- What is the smallest 5 q for which F u (P 0 ,8 q ) is unstable? (|i) 


• Nominal |0. and p-k margins should be similar 

- both methods use the same theoretical plant model 

- both methods use the model with no accounting for uncertainty 


$\mu$ computes a flutter margin for the nominal system by 
considering the smallest perturbation to dynamic pressure which 
causes an instability. Considering all perturbations to dynamic 
pressure allows (I to extrapolate to the flutter boundary from 
a particular stable flight condition. The nominal margins 
computed with the |i and p-k methods should be similar 
since both methods utilize the same theoretical aeroelastic model 
with no accounting for modeling errors. 
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(i Method : Robust Flutter Margins 


• Uncertainty should be included in model 

- Nominal model is not accurate 

- Flutter margins should consider errors in the model 

• Formulate uncertain P 0 at nominal q 0 

- Add feedback operator S q to account for flight conditions 

- Add feedback operator A to account for modeling uncertainty 

• |i computes robust flutter margins 

- Utilize block diagonal uncertainty A={ 8 q ,A} 

- Iterate over scalings on 8 q until p=l 

- (i is worst-case 8 q with respect to A 


The (I method truly becomes a unique valuable tool for flutter 
analysis by accounting for modeling uncertainty. This uncertainty 
is essential since the nominal model is never exactly accurate 
so the flutter margins computed with that nominal model may 
be arbitrarily different than the margins of the true aircraft. 

The basic procedure for including uncertainty is to introduce 
feedback operators to appropriate elements of the system 
and utilize the LFT framework to express an uncertain plant model 
comprised of a nominal plant model and a single structured 
uncertainty set. A process iterating between adjustments to 
norm bounds on the uncertainty and computations of |i can 
be used to compute a flutter margin associated with the 
desired (i=l condition. 
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Consider Pade approximation of Q with lag uncertainty 

s 

y g = Q(s)u q = 5 +|i 0 (l +Ap ) ug 
Separate perturbation from nominal state equation 

X = -$X~J$Uq 

= -pc-O + AfOx-VMl + A p)«, 

= — — -^QoUq + Ap(- |3aX — •^QoUq') 

= - PaX - -J^Uq + ApZ 
= - PaX - ffiollq + W 


Formulate LFT as F u (Q 0 ,Ap) 


X 


-p « 

1 


x_ 


L_ 

r\ 

L_ 


z 

= 

-p„ 

0 


w 



Q o 






0 

1 

- Uq - 

y q 


u q 


An example of modeling uncertainty arises when considering 

aerodynamic lag from a Pade approximation. These lag terms 

may have slight errors so a perturbation can be introduced to 

account for the errors. The perturbation can be separated from 

the nominal dynamics and extracted as an uncertainty operator introduced 

in a feedback manner. 
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(i. Method : Robust Flutter Margins 


• Denote uncertainty A=diag{5 q ,A} 

- Parameterization of Dynamic Pressure 

- Structured Uncertainty 



Nominal Flutter Question : 

- What is the smallest 6 q for which F u (P 0 ,8 q ) is unstable? 

Robust Flutter Question : 

- What is the smallest 8 q for which F U (P 0 ,A) is unstable for some A in A? 


All the individual uncertainty operators, such as the lag uncertainty, 
are combined into a single structured uncertainty operator using 
the LFT framework. The nominal flutter margins are computed by 
considering the smallest destabilizing perturbation to 
dynamic pressure that causes an instability in the nominal plant P. 
The robust flutter margins are computed by considering 
the smallest destabilizing perturbation to dynamic pressure that 
causes an instability in any member of the set of plants F U (P,A). 
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Method : Robust Flutter Margins 


• What does a robust flutter margin really mean? 

- F U (P,A) is a set of plant models 

including nominal and perturbed dynamics 

- Each perturbed plant in the set has a flutter margin 

- Robust flutter margin is the smallest margin 

of all margins for the set of plants 


• Robust flutter margins are conservative 

- Uncertainty can only decrease stability margin 

- Robust |i margin is never greater than nominal p margin 


The robust flutter margin concept seems a significant departure from 

traditional flutter margins; however, the robust flutter margin 

can be interpreted in terms of nominal flutter margins. A unique 

plant model is obtained for each uncertainty operator in the set 

and a nominal flutter margin can be computed for each of these 

unique plant models. The robust flutter margin is simply 

the smallest margin of these individual margins. The operator theory 

and functional analysis concepts are introduced to ensure the 

robust flutter margin analysis considers all uncertainty operators in the set. 
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fi Method : Incorporating Flight Data 


IWSimm- 


• Flight data should be utilized in the analysis 

- Theoretical P may not be accurate 

- Theoretical A may not be reasonable 

• Flight data can be incorporated in several ways 

- Formulate F U (P,A) (difficult) 

- Formulate P (problematic) 

- Formulate A (advantageous) 

• An approach 

- Utilize the original theoretical plant as the nominal model 

- Choose location and structure for uncertainties 

- Use model validation to choose size of uncertainties 


Flight data provides the only true indication of the aircraft properties 
and should be included in the flutter margin analysis. The optimal 
approach would use this data to directly identify parameters in the 
model and an associated uncertainty set but this is extremely difficult. 
Approaches to utilize the data to identify a plant model are valuable 
but problematic since they are unreliable with poor quality data. 

The approach investigated at NASA Dryden Flight Research Center is 
to utilize the flight data to identify the errors and uncertainty 
in a theoretical plant model. 


25 




NASA 


(i Method : Incorporating Flight Data 


• Why formulate only A with flight data? 

- Computational algorithm is straightforward 

• Formulating P via nonlinear optimization is nonunique 

• A can be updated with simple scalar multiplication 

- Method works with poor quality flight data 

• Use a priori theoretical A 0 if data can’t be used for validation 

• A 0 can be used for all or specific frequencies 

- Method accounts for time-varying behavior 

• Difficult to choose best data to formulate P 

• A can be norm bound accounting for all time variations 

- Exact model parameters do not have to be computed 

• Parameter estimation of P is difficult with poor data 

• Model validation is transfer function based 


There are several advantages to using flight data to compute 
the uncertainty description as compared to identifying plant model 
parameters. One main advantage is the ability to utilize data with 
varying levels of noise. Another advantage, which will be utilized 
for the on-line implementation, is the ability to account for 
time-varying dynamics of the airplane. 


26 




|0. Method : Incorporating Flight Data 


Flight data can validate the uncertainty 

so A is a reasonable indicator of modeling errors 

• Choose A large enough to be accurate 

- Increase size of A so no data can invalidate the model 

- Increase size of A at particular frequencies where data differs from P 

- Increase size of particular blocks of A 

• Choose A small enough to reduce conservatism 

- Reduce size of A so some data almost invalidates the model 

- Reduce size of A at particular frequencies 

- Reduce size of particular blocks of A 


An uncertainty description must be chosen that is a reasonable 
indication of modeling errors. This implies a tradeoff must be met 
between choosing A large enough to account for all errors 
but choosing A small enough to reduce conservatism in the 
robust margins. 
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(j. Method : Incorporating Right Data 


• Update robust model through A 

- Finite element model does not change 

- A may be determined from flight data observations 



• (I flutter margins are worst-case 

- [i is worst-case with respect to uncertainty 

- |i is worst-case with respect to errors observed with flight data 


The [i method can be formulated using a theoretical plant model 
and measured flight data. The theoretical plant model is assumed 
to be the best estimate of the aircraft dynamics so this model 
is not changed throughout the analysis. The flight data is 
incorporated entirely through the uncertainty operators. 

(i computes a flutter margin for the theoretical plant 
that is worst-case with respect to the uncertainty and thus 
is worst-case with respect to the observed flight data. 
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|i Method : Robust Flutter Margins 


Include uncertainties to formulate robust model. 

• Uncertainty from Modeling Principles 

- Unmodeled dynamics and nonlinearities 

- Errors in structural elements and parameters 

- Inaccuracies in unsteady dynamic force model 

• Uncertainty from Flight Test Data 

- Measurement of excitation 

- Nonrepeatibility 

• Uncertainty from Signal Processing 

- Waveform basis 

- Assumptions of stationary, linearity, time-invariance 

• Uncertainty from control surfaces on aircraft 

- Nonlinear actuators 

- Hysteresis and Freeplay 


There are many sources of uncertainty in a plant model. The 
obvious sources are errors in the theoretical plant dynamics 
arising from inaccuracies in the structural and aerodynamic models. 
Utilizing flight data to validate the model introduces additional 
uncertainty since the data measurement and processing may not 
be accurate. Also, the aircraft may display behaviors for which 
the measured excitation signal does not account. 
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p. Method : Properties of p 


p has several desirable features as a stability margin 
as compared to parameters such as damping 


- Conservatism is a measure of sensitivity 

• Models sensitive to errors will be conservative to uncertainty 

- (i analysis can determine worst-case perturbation 

• Indicates worst-case flutter mechanism 

• May indicate active and passive flutter control strategies 

- p is a stability predictor 

• Damping is only guaranteed informative at instability 

• \x extrapolates across flight condition 

• fi is linear across dynamic pressures with no A 

• fi is generally well-behaved across dynamic pressure with A 


There are several desirable properties of p that make this 
parameter an advantageous flutter margin as compared to damping, 
p presents more information than damping since the 
conservatism between the robust and nominal margins is an 
indication of modeling sensitivity. Most importantly, p 
is a stability predictor and extrapolates to an unstable flight 
condition whereas damping is merely a stability indicator and 
is only truly informative at the point of instability. 
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(0, Method : Assessment 


• Limitations 

- Model based 

- Linear 

• Difficulties 

- Difficult to choose structure of uncertainty 

- Difficult to choose method of updating the model 

- Difficult to choose flight data used in model validation 


The method is a powerful tool for flutter margin analysis 

but it is only useful when its limitations and difficulties are understood. 

This method is limited to applications of analyzing stability of aircraft 

with a plant model comprised of linear operators. The main 

difficulty lies in choosing the uncertainty description, fi 

is directly related to this uncertainty description so the robust flutter 

margin may be meaningless if the uncertainty description is meaningless. 
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|i Method : Assessment 


• Advantages 

- Can utilize good and bad quality flight data 

- Can use off-the-shelf hardware and software 

- Computational cost is often not extravagant 

- Applied to a real problem and got acceptable results 

• Severed situations where |i method is ideal 

- Want to consider several variations of a model (configurations) 

- Know locations and bounds on model element variations 

- Want some indication of model sensitivity 

- Modeling errors are not huge 

- Want closed-loop margins with a controller (ASE) 


The |i method, despite its limitations and difficulties, can 
be effectively applied to many aircraft applications. This method 
is ideal for considering several variations and configurations 
of an aircraft that only slightly perturb the nominal model. 

This method is also ideal for analyzing aeroservoelastic stability 
since it is a trivial extension to include a controller in the 
formulation and consider the closed-loop dynamics. 
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Worst-Case Flutter Margins : Post-Flight Analysis 


• | i Method can be utilized for post-flight analysis 

- Aircraft demonstrated certain stability properties at test points 

- Flight data should be used to update the model 

- |i computes worst-case flutter margins with respect to A 

• Post-flight analysis with p can be valuable 

- Determine worst-case flight envelope 

- Identify inaccuracies and sensitivities in a model 

- Demonstrate test points that should be repeated 


The p method can be utilized for post flight data analysis 
to compute worst-case flutter margins and determine a 
worst-case flight envelope. The conservatism between the 
robust and nominal margins demonstrates the sensitivity of the model 
and may indicates areas of the flight envelope for which more data 
should be generated due to dramatic differences between the 
theoretical nominal and flight derived robust margins. 


33 



r* 

NA5A 


Worst-Case Flutter Margins : Post-Flight Analysis 


• Available information 

- Nominal aeroelastic model P(M) at different Mach numbers 

- Multiple data sets 

• Multiple test points at identical flight conditions 

• Multiple test points at varying flight conditions 

• Choose uncertainty A from data at certain Mach 

- A(M) varies with Mach number (local) 

- A constant for all Mach numbers (global) 

- A(M) is piecewise-constant (hybrid) 

• Approach : compute |i(F u (P(M),A(M))) for each M 

- Size of A(M) relates model accuracy 

- Conservatism of (A relates model sensitivity 

- Value of Ji relates worst-case flutter margin 


There are several approaches that may be used for deriving an 
uncertainty description from multiple sets of flight data. This 
uncertainty description is formulated as a function of Mach number 
to associate with each plant model describing the 
dynamics at each Mach number. A robust flutter margin is 
computed for the associated plant and uncertainty models at 
each of these Mach numbers and the worst-case flight envelope 
is computed with respect to both dynamic pressure 
and Mach flight conditions. 
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Worst-Case Flutter Margins : A 


• Local approach uses unique A(M) at each Mach 

- Model validation only considers data sets from that Mach 

- Advantage : less conservative 

- Disadvantage : susceptible to poor data at a particular Mach 

• Global approach uses constant A at all Mach 

- Model validation considers data sets from all Mach numbers 

- Advantage : A is worst-case with respect to all test points 

- Disadvantage : conservative at well-modeled test points 

• Hybrid approach uses piecewise-constant A(M) 

- Model validations considers all data from certain Mach numbers 

- Advantage : A is reasonable for conservatism and accuracy 

- Disadvantage : may be difficult to group Mach numbers 


The choice of which data sets to analyze affects the uncertainty 
description. A local approach generates an uncertainty description 
for the plant at a particular Mach number by only considering 
flight data measured at that Mach number. A global approach 
considers all sets of flight data and generates a single uncertainty 
description that is worst-case with respect to variations observed 
throughout the flight envelope. A hybrid approach can be used to group 
sets of flight data together by assuming certain Mach numbers 
at certain test points, such as subsonic flight conditions, are 
strongly related to each other but not to other flight regions. 
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• F/A-18 Systems Research Aircraft 

- well developed finite element model 

- complete set of p-k flutter solutions 


Post flight analysis is used to compute worst-case flutter margins 
for the F/A-18 Systems Research Aircraft (SRA). This aircraft 
is an excellent application for this method since a well-developed 
theoretical model is available along with a complete set 
of nominal p-k flutter margins. 
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Worst-Case Flutter Margins : F/A-18 SRA 



• Flight data from Wingtip Exciters 

- multiple sine sweeps from 3 to 35 Hz 

- symmetric and antisymmetric operation 

• Test points throughout flight envelope 

- Mach range : 0.6 < M < 1 .4 

- altitude range : 10 kft < h < 40 kft 


The F/A-18 SRA is also an excellent application since a large database 
of flight data has been recorded in response to a set of wingtip 
exciters. The test points for these response measurements cover 
a large range of flight conditions throughout the flight envelope. 


37 



Pw A«5A 


Worst-Case Rutter Margins : F/A-18 SRA 


• State-space model aeroelastic dynamics 

- Finite element structural model has 34 modes 

• 6 rigid body modes 

• 14 symmetric elastic modes 

• 14 antisymmetric elastic modes 

- Unsteady aerodynamic model has 84 states 

• 56 states for symmetric forces 

• 28 states for antisymmetric forces 

• Input and Output signals on each side of the aircraft 

- 5 output accelerometers 

- 1 input exciter force 


The state-space model for this system includes 34 modes from 
the structural model and 84 states accounting for the unsteady 
aerodynamic forces. There are 5 accelerometers measurements on 
each wingtip of the aircraft with a measurement of excitation 
force on each wingtip. 
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Worst-Case Flutter Margins 


F/A-18 SRA 


• Uncertainty in Modeling 

- Parametric 

• Mass 

• Damping 

• Aerodynamics 

- Dynamic 

• Nonlinearities 

• Unmodeled dynamics 



Uncertainty in Flight Data 

- Flight Test 

• Exciter Performance 

• Nonrepeatibility 

• Flight Path } 

- Fourier Processing 

• Time-varying transients 

• Nonlinear behavior 




Areas of modeling uncertainty can be seen by visually inspecting 
some flight data sets. The top plot demonstrates two flight data 
transfer functions taken at identical flight conditions. These 
transfer functions show a small variation of approximately .4 Hz 
in the modal natural frequency of the Wing 1st Bending mode. 

The bottom plots demonstrate the concept of nonrepeatibility that 
affects flight flutter testing. These two plots show transfer 
functions generated by response data at identical flight conditions 
for which the modal frequencies and levels displayed are clearly different. 





A a : modal uncertainty 

- 5% in natural frequency 

- 15% in damping 

- 15% in aerodynamic lags 


A m : input uncertainty 

- 10% at low frequency 

- 100% at 150 Hz 

- 500% at high frequency 


An uncertainty description is chosen for the F/A-18 SRA based on 

the observed flight data variations. A parametric modal uncertainty operator 

is chosen to introduce 5% uncertainty in natural frequencies, 

15% uncertainty in damping ratios and 15% uncertainty in 
aerodynamic lags. A complex dynamic input multiplicative uncertainty 
is included to introduce 10% uncertainty at low frequencies and 
increasing to demonstrate the model is poor at high frequencies. 

Sensor noise is included along with the perturbation to 
dynamic pressure that allows JlI to extrapolate to the flutter margin. 
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Worst-Case Flutter Margins 


F/A-18 SRA 



Nominal Symmetric Flutter Points 


The initial calculation of flutter margins considers only the 
nominal plant dynamics with no consideration of the uncertainty 
description. This plot shows the p-k margins as solid 
lines and the nominal p margins as circles. These values 
are quite similar to indicate the p method can accurately 
compute flutter margins. The frequencies of each unstable flutter 
mode are not given on this plot but are shown in the published reference 
documents to match closely. 
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Worst-Case Flutter Margins : F/A-18 SRA 



Symmetric Flutter Margins 


Robust flutter margins are computed for the symmetric modes 
using p with respect to the uncertainty description. 

The dashed line indicates these worst-case margins which are 
more conservative and closer to the flight envelope than the nominal 
margins. The robust flutter margins at transonic flight conditions 
are particularly conservative. 
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Worst-Case Flutter Margins 


F/A-18 SRA 



Antisymmetric Flutter Margins 


Robust flutter margins are also computed for the antisymmetric 
modes and are shown by the dashed line in this plot. These 
margins show a similarity to the symmetric mode margins in 
that the transonic flight condition is especially sensitive 
to modeling uncertainty. 






Worst-Case Flutter Margins : F/A-18 SRA 


• Features of Symmetric 

- Similar margins at supersonic 

- Worst-Case at M=1 .2 close to 15% margin 

• Features of Antisymmetric 

- High sensitivity at M=1.6 

- Subcritical hump mode at M=.95 


• Features of Both 

- Poor nominal match at M= 1 . 1 

- High sensitivity at M=1 . 1 

- Worst-case margins outside 15% margin 




The entire set of robust flutter margins can be considered by 
evaluating the margins for the symmetric and antisymmetric modes. 

These margins, despite the conservatism associated with including 
uncertainty, are more than 15% in airspeed from the flight envelope of the 
F/A-18 SRA. 
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Flutterometer : On-Line Implementation 


• Extend method to on-line predictions during flight 

- Update uncertainty description at each test point 

- Update flutter margin at each test point 

• Flutterometer displays on-line predictions 

- p prediction directly accounts for flight data 

- p prediction extrapolates to flutter boundary 

• Flutterometer can improve flight test efficiency 

- Test point data tracks time-varying dynamics 

- Worst-case margins provide confidence 

- p is a better behaved stability margin as compared to damping 

- Confident margins can be used to adapt flight plan 


The (I method for flutter analysis is extended to 
an on-line implementation by considering a test point approach. 
Flight data is gathered at each test point and immediately 
used to generate an uncertainty description and an associated 
robust flutter margin. The robust flutter margin information at 
each test point is displayed via the flutterometer tool which indicates 
the distance to flutter for that test point. This approach allows the 
flutterometer to track time-varying dynamics of the airplane since 
the uncertainty description is continually updated as flight 
data is measured and the p. margin accounts for the time-varying 
dynamics through this time-varying uncertainty. 
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Flutterometer : On-Line Implementation 


NA^A 


• Several methods to choose A during a flight test 

- A considers data from current test point (local) 

- A considers data from all test points (global) 

- A considers data from recent test points (hybrid) 


• Approach is similar to current methods 

- Trim at a stabilized test point 

- Record transfer function data 

- Transfer data to analysis computer 


< Estimate modal damping 

Determine reasonable uncertainty A 
Compute worst-case flutter margin p 


- Determine conditions for next, if any, test point 


The (I method can be implemented in a manner similar to computing 
damping estimates. The flight data can be simultaneously analyzed 
by traditional methods and the jli method on separate computers. 

The flight conditions for the next test point can be determined 
using the combination of information from the traditional damping 
estimate and the new j x method. 
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Flutterometer : On-Line Implementation 


• Local approach uses current data 

- Only data from current test point is used 

- Advantage : less conservative 

- Disadvantage : susceptible to poor data set 

• Global approach uses all data from the flight 

- Model validation considers data from all previous test points 

- Advantage : worst-case with respect to range of flight conditions 

- Disadvantage : may be overly conservative 

• Hybrid approach uses recent flight data 

- Model validation considers current and recent data 

- Advantage : reasonable for conservatism and accuracy 

- Disadvantage : may be difficult to define recent 


The choice of flight data to utilize at each test point allows 
flutter margins to be computed which are worst-case with respect 
to different uncertainty operators. A local approach uses only 
data from the current test point while a global approach uses all 
data from previous test points. A hybrid approach uses a forgetting 
factor and only uses data from recent test points. 
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The flutterometer based on the |i method can be easily 
implemented using current flight flutter test procedures. The 
only change is an additional analysis operation using the recorded 
flight data that computes a robust flutter margin. The resulting 
decision on future test points is considered confident since 
the damping estimates are accurate indicators of the current stability 
properties and the p method provides significant additional 
information about the distance to a flutter instability. 
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Flutterometer : F/A-18 simulation 


• Simulate flight flutter test of F/A-18 SRA 

- constant Mach dive at M= 1 .2 

- test points every 100 lb/ft 2 of q 

• Test point procedure 

- Record frequency response data 

- generate uncertainty description 

• Validate current uncertainty levels 

• Increase if necessary 

- Compute worst-case p flutter margin 

• Time is variable in the simulation 

- computation time to validate A 

- computation time to compute p 



A simulated flight flutter test of the F/A-18 SRA demonstrates 
the flutterometer concept. The aircraft is undergoing a constant 
Mach dive at M=1.2 to expand the envelope. Test points are chosen 
at every 100 units of dynamic pressure to illustrate the on-line 
computations with some detail. The time spent at each of these 
test points is determined entirely by the amount of computation 
time required to compute an uncertainty description and an 
associated robust flutter margin with |i. This plot shows 
the value of dynamic pressure throughout the simulation with 
the length of the horizontal lines indicating the time spent 
at a test point for which the flight condition did not change 
while |! is being computed. 
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Flutterometer : F/A-18 simulation 


• Plant model is not accurate 


- error in structural damping 

- error in initial mass value 

- unmodeled time-varying mass 


(true aircraft is 10% higher) 

(true aircraft is 95% of heavyweight) 
(true aircraft varies 5% in 20 minutes 


• Modal uncertainty increases as plant dynamics change 


- Modal parameters are time-varying 

- A increases as mass decreases 



A time-varying theoretical plant model of the F/A-18 SRA is used as the 
true plant model while a variation of this model is used as the 
nominal dynamical model. The nominal model has an error in structural 
damping and does not account for the time- varying mass of the true 
aircraft. The plot shows the increase in modal damping uncertainty that 
is required throughout the simulation to ensure the flight data 
recorded at each test point does not invalidate the uncertain model. 
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Flutterometer : F/A-18 simulation 



The top plot shows the true and computed flutter margins 
for each test point. The dashed line is the flutter margin of 
the true aircraft and decreases with time due to the time-varying 
mass of the aircraft. The solid line displays the robust flutter 
margin and also demonstrates a time-varying behavior due to the 
calculation of an uncertainty description and |i at each 
test point. This robust flutter margin utilizes uncertainty 
to account for the time-varying mass and remains conservative 
to the true flutter margin throughout the simulation. The dotted 
line near the top of the plot is the nominal flutter margin which 
does not account for any uncertainty in the nominal model. 

The bottom plot shows the modal damping ratio at each test point 
for the mode that goes unstable at the M=1.2 flight condition. 

The damping remains fairly constant until minute 17 when it 
decreases sharply to indicate the oncoming instability. 
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Flutterometer : F/A- 18 simulation 


•M = 1.2 

Consider minute 17 test point : *<1 = 1700 

•h = 6147 ft 



Damping trend indicates potential impending instability 
but flutterometer quantifies distance to instability 


These plots demonstrate the information displayed in the control 
room at the minute 17 test point. The flutterometer on the left 
indicates the aircraft can drop 7500 feet before a flutter condition 
is encountered. The damping trend on the right indicates an 
instability may be near but the flight conditions associated with 
that instability can be not accurately determined from this trend. 
These plots clearly show the usefulness of a valid flutterometer 
tool for flight flutter testing. 





Flutterometer : F/A-18 simulation 




Computations used standard equipment 


- Pentium 200 MHz computer 

- Matlab and p-Tools software 

CPU time was reasonable 

- time to validate A 

- time to compute ji 



What affects CPU time? 


simulation time (min) 


- Number of uncertainties 

- Number of validation iterations 

- Number of frequency points to compute p 

- Number of states in the model 


This simulation uses standard hardware and software often 
available in a flight data analysis facility. The computational 
time required at each test point is reasonable and does not 
introduce an excessive burden on a flight test program. The plot 
uses a star symbol to indicate the CPU time required 
at any test was not greater than 180 seconds. 






Flutterometer : Assessment 


Flutterometer has same limitations/advantages as |X method 
with additional issues due to the implementation 


• Limitations 

- Needs stabilized test points 

- Needs constant Mach envelope expansion 

- Needs transfer function data (output and input) 

- Needs some computation time (not real -time) 

• Advantages 

- Tracks time- varying dynamics and flutter margins 

- Computational time is reasonable 

- p is better behaved stability parameter than damping 


The flutterometer tool is based on the p method and 
has the same limitations and advantages as the p method with 
several additional considerations. This tool is only directly 
applicable for flight programs using stabilized test points that can 
vary dynamic pressure while keeping Mach constant. The flutterometer 
can be extremely useful if these limitations are not problematic 
since the time-varying flutter margins can be tracked and 
p is a much better behaved stability parameter than damping. 
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Conclusions 


• [X method is improvement to traditional methods 

- Flight data and model are utilized 

- Worst-case flutter conditions are determined 

• Flight data can be utilized in different ways 

- Post-flight analysis to update stability margins 

- In-flight analysis to update stability margins 

• Flutterometer can improve flight test efficiency 

- Test point data tracks time-varying dynamics 

- Worst-case margins provide some level of confidence 

- p replaces poorly behaved damping as a stability margin 


The p method is a significant improvement to traditional 
flutter margin analysis methods since it directly accounts 
for modeling uncertainty by including flight data. This method 
can be used in several ways including an on-line implementation 
to develop a flutterometer flight test tool. This flutterometer 
can improve flight test efficiency by providing a confident 
measure of the distance to flutter so test points can be chosen 
that quickly and safely expand the flight envelope. 
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(X Method : Research Extensions 


• Aeroservoelastic stability margins 

- flutter and buffet stability margins 

- gain and phase controller margins for elastic plant 

• Reducing conservatism in stability margins 

- Wavelet filtering to reduce conservatism 

- Model updating to reduce conservatism 

• Analyzing nonlinear and LCO dynamics 

- Using |i to predict certain classes of transonic LCO 

- Extending LFT framework to include nonlinear operators 

• Demonstrations and Applications 

- Experimental Flight Testbeds 

- Flight Test Programs 


There are several research extensions to the |X method that 
are under investigation. The computation of aeroservoelastic 
stability margins is performed for an F/A-18 aircraft and results 
will be published in the 1998 AIAA SDM conference. Methods of 
reducing the conservatism in the robust margins are being considered 
using techniques such as wavelet filtering and model updating. 

The fl method is being extended to account for nonlinear 
limit cycle oscillations by including nonlinear operators in 
the LFT framework. Also, the (I method is being applied 
to several flight test programs to validate and improve 
the implementation issues. 
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The ^ method is being applied to several small order testbed 
systems. These systems can be flown at the flutter boundary 
to ensure the flutterometer predicts the unstable flight conditions. 
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The |i method is also being considered for several aircraft 
flight flutter test programs. The cost of these programs could 
be dramatically reduced with even a small increase in flight 
test efficiency using the flutterometer. 
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|i Method : Discussion Topics (general) 


• The (I method 

- Why is p better than traditional/new methods? 

- What do robustness and p actually mean? 

- How does p relate to damping? 

- How do I develop an uncertainty description? 

- What level of confidence can I place in these margins? 

- Does p have any usefulness in analyzing nonlinear systems? 

• Incorporating flight data 

- What types of data can be used? 

- What is the model validation actually doing? 

- What are local and global updating schemes? 

• Flutterometer 

- What flight test and control room procedures can be used? 

- Can this be done in real-time? 

- What if flutterometer reads “0” margin at my current test point? 




p Method : Discussion Topics (Boeing) 


• Control room issues 

- Hardware and software requirements 

- Interaction with other Boeing tools and procedures 

- Work level of the user during a flight test 

• Utilizing flutterometer 

- Availability of algorithms 

- Cost and effort to develop and implement 

- External funding possibilities and interests 
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